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ON FLUCTUATIONS OF RIEMANN'S ZETA ZEROS 
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O ' Abstract 

(N 

«i . It is shown that the normalized fluctuations of Riemann's zeta zeros around 

fn ' 

their predicted locations follow the Gaussian law. It is also shown that fluc- 
tuations of two zeros, 7^ and jk+x, with x ~ q (log ky , /3 > 0, for large k 
follow the two-variate Gaussian distribution with correlation (1 — /3)_,_ . 

1. Introduction 

This paper is concerned with the statistical properties of the Riemann zeta function ze- 
ros. This subject originated in 1944, when Selberg [15 1 showed that the number of zeros 
in a sufficiently long interval on the critical line can be described by the Gaussian law (see 
also |fT6l . ||5l , ||9^ , |[8l). In the 1970s, Montgomery and Dyson discovered the remarkable 
fact that the spacings between the zeta zeros resemble the spacings between the eigenvalues 
^ ■ of random Hermitian Gaussian matrices. This resemblance was substantiated analytically by 

Q ' Montgomery ifTTTl and supported numerically by Odlyzko |[T3l (see also lfT4l . f2\). The con 

cn ■ 
vn 

en 

O ' tion by showing that at the mesoscopic level Riemann's zeros have correlations previously 

found by Diaconis and Evans |i4J for eigenvalues of unitary random matrices. 

The motivation for this study comes from a recent paper by Gustavsson Q, who showed 
\Ji • that eigenvalues of random Hermitian matrices fluctuate according to the Gaussian law. Our 

^ . goal is to investigate the statistical fluctuations of Riemann's zeros around their predicted 

- - - positions and to show that these fluctuations also follow the Gaussian law. 

We assume Riemann's hypothesis and denote the non-trivial zeros of Riemann's zeta func- 
tion by 1/2 + ijk- We consider only zeros with positive imaginary part, 7^ > 0, and order 
them so that the imaginary part increases, 71 < 72 < • • • • 

Let J\f{T) denote the number of zeros with the imaginary part strictly between and T. If 
there is a zero with imaginary part equal to T, then we count this zero as 1/2. 
Define 

5(r):=ilmlogC(J+ir), 



nection between zeta zeros and random matrix eigenvalues drew much attention, as can be 
seen for example from review papers in fTOl. Recently, Bourgade [Ij supported this connec- 
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2 V. KARGIN 

where the logarithm is calculated by continuous variation along the contour a + iT, with a 
changing from +00 to 1/2. 

It is known (see Ch. 15 in fSl) that 

M{T) = ^\og^ + l + S{T) + 0^ ^ 



27r ° 27re 8 ^ ' \l + T 

Let tfc be the solution of the equation 

t t 7 , , 

— log h- = k- 1/2, 

27r ^ 27re 8 ' ' 

and let 

V21oglogtfc 

Let 7fc denote the imaginary part of the A;-th Riemann zero and define 

h = ^-^^^^. (1) 

The quantities /^ measure fluctuations of Riemann's zeros from their predicted locations 
tfc. In order to study the statistical properties of /^ we introduce a probability space {Q.,B,¥}, 
where Q. = [Q,\] ,B is the u-algebra of Borel subsets of $7, and P is the Lebesgue measure on 
B. 

Let us fix a > and 61 G (1/2, 1]. We will denote [aN] by Kn, [N'^\ by Hn and the 
interval of integers [Kj^, Kn + H^ — 1] by In. We define a sequence of random variables 
fW by the following formula: 

/W(a;):=A,(^,,), (2) 

where 

k{N,u;):= [Kn+u;Hn\. 

( [xj is the largest integer which is less than or equal to x.) Hence k {N, to) is a random variable 
taking its values in integers and uniformly distributed on the interval In- Note that 

p|/(^) G (a,6)| = -^\{k:kG lN,fk G {a,b)}\ , 
L J Hn 

and 



Hn 



We will prove the following theorem. 

Theorem 1.1. Suppose that random variables f^ ' are defined as in (O with parameters 
a > and 1/2 < 6 < 1. Then, for every real S^, as N ^ 00, 



{'"">^)-^Jr-''"'- 



The following simple Corollary of Theorem II. H is also useful: 
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Corollary 1.2. For every p > 0, 



lim E f/(^)y = -^ r xPe--'/'dx. 



It is also interesting to ask how fk and fy are related when k and k' are sufficiently close 

fine 

/r^H=4(^,..), ^=1,2, (3) 



to each other. More precisely, define random variables /{ and /2 by the formula 



where 



ki {N, u) = N+ [ujN\ , 
k2 {N, io) = N+ [ujN\ + \a (log N)^ 



and a > 0, /3 > 0. 



Theorem 1.3. Suppose that random variables f^ are defined as in with parameters 
a > and /? > 0, and suppose that Xi, X2 are zero-mean Gaussian random variables with 

E [Xf) =land¥. (X1X2) = (1 - /3)+ . Then as N -^ 00, 

(i) the joint cumulative distribution function of f /} , /2 j converges pointwise to the joint 



cumulative distribution function of {Xi,X2), and 

(ii) the joint moments 0/ ( /{ , /2 ) converge to the corresponding joint moments of [Xi , X2 ) 

The rest of the paper is organized as follows. Section |2] outlines the scheme of the proof 
of Theorems 11.11 and 11.31 Section [3] introduces some technical tools that we will need in 
the proof of the main theorems. Section |4] proves a modification of the key approximation 
result by Selberg. Section |5] calculates the moments of the approximate function Sx- Section 
[6] calculates the moments of S and concludes the proof of Theorem ll.il Section |7] proves 
Theorem 12.31 Section [8] proves Theorem 124] and concludes the proof of Theorem 11.31 And 
Section |9] concludes . 



2. Outline of Proofs 
It is convenient to define 

V27rS{tk + ^ak) 

A^ • — ^ ) 

Vlog log tk 
and a corresponding sequence of random variables 

X(^)M:=Xfc(;v,.). (4) 
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A connection between X^^' and f^^^ can be seen as follows. For every real ^, 

\{k : k e In, M (tk + ^ak) < k - l/2}\ 



Hn 
1 

Hn 
1 

Hn 

1 

Hn 

1 



k:kelN,^-^^^^\og^-^4^ + l + S{tk + i^k) + 0{l/tk)<k-l/2 



27r 



27re 



^ : ^ E /., 5 (t, + C.,) < -i,r^^^^ ( 1 + ^ K « (1 A.; 



27r2 V log^A: 

log27r 



k:k£lN, ,, ■ ^ — < -n 1 + - — — + o (1/tfcj 



Since f ^ is asymptotically close to 2irk/ log A;, it follows that 

(5) 
Hence for large N, the distribution of the random variable /^^^ is essentially determined by 
the distribution of the random variable X^^\ 

In this connection, it is appropriate to recall the following theorem by Selberg (Theorem 3 
in lUl). Let 

V log log t 

Theorem 2.1 (Selberg). Assume RH, and let T'' < H < T'^, where a > 0. Then for every 
k> 1 

2k\ 



+ 



logN 

N 



1 pT+H 

— / \X {t)\'^'' dt = —j: 

H Jt k\2k 



+ 0(1/ log log T), 



with the constant in the remainder term that depends only on k and a. 

In other words, the even moments of the function X (t) behave as the moments of a stan- 
dard Gaussian variable. This was refined in [5 1, where it was shown in particular that for every 
interval /, 

1 £^'' 1, [X (t)] dt = ^l^ e-^'/'dx + 0(1), 

where 1/ denotes the indicator function of interval I. We prove a modified version of this 
result. 

Theorem 2.2. Suppose that random variables X^^' are defined as in (]?]) with parameters 
a > and 1/2 < ^ < 1. Then, for every real ^, as N ^ oo, 



El(_,_^] (X(^) 



1 



-^'/'dx. 
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Theorem 11.11 immediately follows from Theorem I2.2l and inequalities ^. We will prove 
Theorem l2.2l by the method of moments, which says that in order to establish the convergence 
of a sequence of r.v. in distribution to the Gaussian law it is enough to show the convergence 
of every moment (Example 2.23 on p. 18 in van der Vaart IITtI ). That is, it is enough to show 
that 



E f X(^) 



1 



X e 



'^"'^dx 



(6) 



V^ J- 

for every integer r > 0. We will show this in Section[6]in Corollary [67 

Corollary 11.21 follows because © and ^ imply that [f^^') are asymptotically uni- 
formly integrable for every n > and therefore the moments of /^^^ converge to the moments 
of the limiting Gaussian distribution (see Theorem 2.20 in van der Vaart [17l). 
In order to prove Theorem 11.31 define random variables 



.(TV) 



Xl'"{uj):=Xk^^^^^),i = l,2, 



(7) 



where ki {N, uj) are as in Q. Let us use notations 

l[^'> := [N,2N - 1], 4^^ := \n + [a(logiV)'^ 
(si,S2) . Then, 

AN) , . AN) ^\ 



, 2iV - 1 + 



Q(logiV)^jj , and s 



iV" 



N- 



N- 



{s : S, E 4''\M{ts,+C^(7sJ} <h- 1/2 
(N) V^TtS (tsi + ii(T, 



That is. 



/r>6,/f^>6}< 



s : Si G /; 



xr < -C. ( 1 



i/log log ts 



C\ 



logiV 



and 



'{/f)>ei,/f^>6}>p{xr^ 



< -ii 1 + 



C2 



loffiV 



+ 



+ 



< -Ci u + 



logiV 



N 



logiV 

N 



log27r 
log ts. 



(8) 



(9) 



In words, the joint cumulative distribution function of /{ and /g approaches that of 
Xf^andxf). 

First of all, we have the following result for the random variables X\ and X^ 

Theorem 2.3. Let X- ' be defined as in ([71). Then, 

:i-/3)+:= 
More generally, the following result holds. 



lim EXf )xf ) 



1 - /?, /f /3 G (0, 1) , 
0, ifii > 1. 



.{N) 



AN) 



.{N) 



Theorem 2.4. Let Xl ' be defined as in (El). Then for every /, m > 0, E \X\ ' \ f X 

converges to E {Xi) {X2)^ where (Xi, X2) is a zero-mean Gaussian random variable with 
E {Xf) =land¥. (X1X2) = (1 - /3)+ . 



+ o{l/ts 
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Theorem 12.31 is a particular case of Theorem 12.41 However, we will prove it separately, 
since its proof is more transparent and shows how the proof of the more general Theorem 12.4 
proceeds. 

Given Theorem |2.4[ we can prove Theorem 11.31 

Proof of Theorem ll.3t Theorem 12.41 implies that the cumulative distribution function 
of ( X| ,^2 j converges pointwise to the cumulative distribution function of (Xi,X2). 
The first claim of the theorem follows immediately from this fact and inequalities ^ and 
(|9l). In addition, Theorem 12.41 and inequalities dS) and Q imply that for all integer a,b > 

0, the random variables ( /} J ( /2 j are asymptotically uniformly integrable. Hence, 

their expectations converge to the corresponding expectation of the limit, E (Xi)"' {X2) (by 
Theorem 2.20 in van der Vaart IITtI ). This completes the proof of the second claim of the 
theorem. D 

The proof of the convergence of moments of X^^') follows the plan of the argument in 
Selberg ifBl . 

Recall that X^^^ is a rescaled version of S (gk) ■ The first step in Selberg's proof is to show 
that S (t) can be approximated by Sx (t) , where 



sin (tlogp) 

Y]- ^ — / 

That is, Selberg shows that 



p<x'^ 



1 rK^H 

-j^ {S it) - SAt))''' dt 

is small provided that K and H are sufficiently large and that x ~ K^ with a sufficiently small 
e > 0. In our case we will need to modify this result in order to show that the integral can be 
replaced by a sum over a discrete set of points. 

The next step in Selberg's proof is to calculate the moments 

fK+H 

H 



1 pK+H 



Again it will be necessary to prove a corresponding result for a sum over a discrete set of 
points. 

Given the results in these two steps, it is relatively easy to calculate the moments of S [t] . 
This will be done essentially as in Selberg's paper. 

3. Exponential Sums 

The changes in Selberg's proof make it necessary to estimate certain exponential sums. 
The main additional tool that we use to handle these sums is the following theorem by van der 
Corput (Theorem 2.2 in [6|). Let e (/ (n)) denote exp [27rif (n)] . 
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Theorem 3.1 (van der Corput). Suppose that f is a real valued function with two continuous 
derivatives on interval I. Suppose also that there is some A > and some n>l such that 

A< \f"{x)\ < k\ 

on I. Then, 

In order to apply this theorem in our situation, we need to estimate derivatives of a function 
g{x) that we are about to define. Let t {x) be the functional inverse of the function 

on interval [to;C«) where to is sufficiently large. Note that t{x) is an increasing concave 
function. Let 

g {x) := t [x) + C — — , (11) 

logt [x) 

where .^ is a real constant. This function is well defined for x greater than some numeric 

constant xq. For x between and xq, we define g{x) in an arbitrary fashion such that g{x) has 

a continuous 3-rd derivative for all x > 0. 



Lemma 3.2. When a; — )• oo, 



and 



11/ ^ 27r 
9 (^) 7, 72 

x (log x) 
27r 



g'" (x) 



x"^ (log x) 



2- 



Proof: By differentiating 



we obtain: 



t(x)^ tix) 7 



27r = t' log (t/27r) , (12) 

," = ^0^ (13) 

t(logt/27r)' ^^ 



and 



tV(-3-log(t/27r)) 
t(logt/27r) ■ 



(14) 



From (ITOl ) and (IT2l ). we get 



X 1 

t ~ 27r: and t' ~ 27r- 



log X log X ' 

and then ([T3]l and ([Hll imply that 

„ 27r , ,,, 27r 

t" o andt" 



2 """ " o n ^2■ 



X (logx) x^ (log x 
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If 



h:-- 



V^log log t (x) 



logt (x) 
then a calculation shows that h" = o{t"), h'" = o{t"') and therefore 



n 
9 ~ 



27r 



X (logx)" 



and g" 



2tt 



x^ (log X 



,2' 



D 



In the following we will use notation g^ for g (k) 
Lemma 3.3. Let 



= log 



Pl+l ■■■P2n 



Pi ■■■Pi 

where 1 < / < 2ri, {pi, ■ ■ ■ iPi} / {pi+i-, ■ ■ ■ P2n} and primes pi < y for all i. Assume 
1< H <cK. Then, 

k=K 



O Hn- 



Kl/2 log K 

Proof: From the assumption, we obtain 

c/y" < \9\ <2nlogy. 
(In order to see the first inequality, let / < n. Then 

1 



+ y"/^K'/HogK). 



Pl+l^ 


■■P2n 





Pi- 


■■Pi- Pl+l ■ 


■P2n 


^ ' 


Pi- 


■Pi 






P1---PI 




yn 



by the uniqueness of integer factorization, and the desired inequality follows. The case I > n 
is similar.) 

Hence, by using Lemma [T2l we find that 



A< \Og"{x)\ < kX 



with 



and 



y'^K (log K) 



2' 



K = 0(ny"logy) 
By applying van der Corput's theorem, we obtain 



K+H-l 



"/2 log 



y e^'a, = o Hn^L^^^JL + yn/2Kl/2 1 ^ 



k=K 



D 
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Lemma 3.4. Suppose 1 < ciK^ < H < C2K, where 9 > 1/2 and ci,C2 > 0. Let y < 

261-1 , 

K 3r ^ ^ and assume that 

I I . log p „ 

\ap\ < A- forp < y. 

logy 



Then, for every integer r > 0,we have 

K+H-l 



E 

k=K 



v<y 



Or 



pl/2+igt, 



2r 



Proof: We can write 



where (3n < A^. Hence, 



0{H). 



Pn 



K+H~l 

E 

k=K 



Es 



ttr 



p<y 



pl/2+igk 



2r 



Z^ pl/2+igk j Z^ j^l/2+igk 

\p<y ' / ra<y'' 



a -5 K+H-1 



J van ^-^ V n 

m,n<y'~ * k=K 



< 



1/3^ 



ify^^ + 2 y 



n<y^ 



m<n<y^ 



\PmPn\ 

Jran 



K+H-l 

E (j^r 



k=K 



The first sum can be estimated as follows: 



■y \Pn\ ^ 



1/3. 



n<y^ 



n 



-^'E^s-^' E^ =0(1) 



n<y^ 



\p<y 



where we used Mertens' result ^p<„ -^^^ = 0(log y) in the last step. 
In order to estimate the second sum we note that 

1/y^ < log |n/r7T,| < rlogy; 

hence we can apply van der Corput's theorem and estimate 



K+H-l 

E 

k=K 



■m\^9k 

n 



<0\ Hr 



y''!'^ \o% 



K^/"^ log K 



y +y-/2^1/2log^ . 



Besides, 



V^ \PmPn\ , 



Iran 



'm<n<y'~ 

By assumptions about H and y, it follows that 

K+H-l 



i^M'-''''' 



E 



\PmPr, 



m<n<y^ 



/mn 



E 

k=K 



m\Wk 
n 



f/"^ loE 



0{H). 



D 
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Lemma 3.5. Suppose 1 < ciK^ < H < C2K, where 9 > 1/2 and ci,C2 > 0. Let y < 

\ap\ < Aforp < y. 



K 3r ^j and assume that 
Then, 



K+H-l 

E 

k=K 



p<y 



Otv 



P 



l+ifffe 



2r 



OiH). 



The proof of this lemma is similar to the proof of the previous one. 



4. An analogue of Selberg's approximation formula 



Recall that 



5.(0 :=—E 



p<x'^ 



sin (tlogp) 



Our immediate goal will be to prove the following result. 

Proposition 4.1. Suppose 1 < ciK^ < H < C2K, where 6 > 1/2 and ci,C2 > 0. Let 
X < K 6n ^. Then, we have 



K+H-l 



E \S{gk)- S.j,{g, 



k)\''' = 0(H). 



k=K 



Proof: Let A (n) = log p, if n is a power of the prime number p, and A (n) = 0, otherwise. 
Also, define 

A (n) , for 1 < n < x, 

^-^^'^ iorx<n<x\ 



Theorem 1 in 1 15] states that if t > 2, 4 < x < t^, and ai = 1/2 + 1/ log x, then 



s-ft) = _i y^ ^kM!i!ilil^l!S!^ + o ' ^ 



vr ^ — ' n°"i logn 



logx 






n<x^ 



+ 



logx 



This gives the following estimate: 
S{t)-S^{t) = o 



y. A(p)-A.(p)pV2--^ ^_,, 



+0 



+0 



+0 



E 

p2<x2 

1 

logx 

log^ 
logx 



V^logp 
A. ip^) 



p2(ai+it) iQgp 



+ 



E 



Ax (P) 



p<x-' 



E 



p 



,2(ai+it) 



logx 



p 



,cri+it 



\r>2 p 



r Ip^ 
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It follows that 

K+H-l 



' K+H-1 



E \S{9k) - SAgk)\''' = 0[ Y. 



k=K 



k=K 

(K+H-l 
E 
k=K 

(K+H-l 
E 
k=K 

'K+H-l 

+ o\ E 

k=K 



y A(p)-A.(py/^-"i ^_,^^. 



p<x'' 



E 

p2<x^ 



V^logp 



^-(P') -,-2.9. 



2n.N 



p2o-i logj5 



P 



2n\ 



V ^/^y p-^'?^ 



p<x 



' (logx)p'' 



2n\ 



E 



Axb') ^-2ig, 



(log x)p2°"i ' 



2n\ 



+ OiH). 



By applying Lemmas 13.41 and 13.51 we find that each of the sums on the right-hand side is 
O {H) . For example, for the first term we can apply Lemma [33] with 

A{p) - A^{p)py^-'^^ 

Op = . 

logp 

The hypothesis of the lemma holds because for p < x we have 



Or. 



I _ g-logp/logx 



<A 



logx' 



and for X < p < x^ the estimate is similar. 
This proves that 

K+H-l 



E l^idk)- Sx{g. 



'kn''' = 0{H) 



k=K 



D 



5. Moments of \Sx {t)\ 
In the next step we show that the moments of Sx {t) are approximately Gaussian. 

Lemma 5.1. Suppose 1 < ciK'^ < H < C2K, where 9 > 1/2 and ci,C2 > 0. Let x < 

29-1 _ 

K dn ^. Then, for every integer n > 1, 
(i) 



K+H-l 

k=K 



and (ii) 



I'" = T^rS^ (^ (i°g i°g ^)" + ^ (^ (i°g i°g ^)""' 

(2tt) n! ^ ^ 



K+H-l 

Yl Sx{gkf''-' = OiH). 

k=K 
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Proof: First, we can write 

where 

V = ri{t)= Y^p-"^-''. (15) 

Hence, 

K+H-l 2n . s K+H-l 

k=K ^"^^) 1=0 ^ ^ k=K 

Here 

K+H-l , K+H-l / ^ ig^ 



k=K p,<x2 V-Pl---^2n ^^^ V 1^1... M / 

If {pi, . . . ,pi} / {pz+1, • • • ,P2n} , then by using Lemma [331 we obtain 



Since 



2n\ 



S;^fis-( S;;.' '-'^^"^ 



hence in the case / 7^ n, we have 

/r+H-i 



k=K 

HI = n, then 



E'^(^^)'^(^^)"" = o(-"(^-^^+-"^^''i°g^) 



o(i/). 



E ^(5.r^(5fcr = H Y. -^—+o{h) 
tK t:^ p^---p- 



Pl...Pn=Pn + l---P2n 



n 



Hiy-] +oin\Hy i \+OiH), 



where the second equality follows from the fact that the number of ways in which a number 
of the form pi . . .pn can be written as a product of n primes is equal to n! if the primes are all 
different and less than n\ if two or more of the primes are equal. Hence, 

K+H-l 
E v{9kTvi9kT = n\H (loglogx)'^ + 0(n\H (logiogx) 



\n—l 
'n\9k) 'I {9k} = n:ii [^Lugiugx) -r <^ yir-n {iogiogx 

k=K 

= n 



ilH (log log Kf + O (nlH (log log E:)"^^) . 
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It follows that 

k=w v^^/ ;=n ^ ^ h=K 



k=K ^ ' 1=0 ^ ' k=K 



(2n)! 



H (log log Kf + 0(H (log log K) 



(27r)2"n! 
The proof of (ii) is similar, except that in this case it is always true that {pi, . . . ,pi} / 

{pi+l,...,P2n-l}-0 

6. Moments of \S {t)\ and conclusion of the proof of the main theorem 

Theorem 6.1. Suppose 1 < ciK^ < H < C2K, where 9 > 1/2 and ci,C2 > 0. Then, for 

every integer n > 1, 
(i) 

E \S idk)?'' = 7^^, (i7 (log log i^r + 0(i7 (log log Kr-1/2)) , 

and (ii) 

K+H-l 

\n-l 



j; {S{gk)f''-' = 0(H{\oglogK) 



k=K 

26-1 



Proof: Take x = i^"%^"^. Then, 

l-S" {9k)f"' = \Sx {9k) + {S {gk) - Sx (Sfc))!^" 



\Sxi9k)?'"' + I ^\Sx{9k)\^'^ ''\S{gk) - Sx{9k)\'' I 



By using the Holder inequality, we can estimate 

K+H-l /K+H-l \ l-'^/2n. /K+H-l \ ^Z^" 

E i^-(5fc)i'"-n5(5fc)-5.(5fc)r < E i^-(5fc)i'" E 1^(5/.) -^.(5^)1'" 

A;=K V fc=ii' / V k=K J 

= 0([H (log log K)«] 1-^^/2" i7^/2n] =0(h (log log K)"-l/2 



2n _ K^ny. / „ _Ln fn nr.<r ir.<r A-^"- 1/2 



Hence, 

E \S{9k)f'' = -^^(H{loglogKr + 0(H{loglogK) 
For the proof of (ii), we estimate 

(2n-l \ 

Y.\SA9k)\''''''-''\S{gk)-SA9kWy 

and note that 

2n—l-v 1+v 

K+H-l /K+H-l \ 2^ /K+H-l \ ~2^ 

E \SA9kT"'~''\S{gk)-SAgkW<[ E l^-(ff't)l'" E l'9(5-fc)-'5x(5fc)l^'" 

k=K V fc=_ft: / V k=K J 
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where we used the Holder inequaUty with p = 2n/ (2n — 1 — u) and q = 2n/ {1 + v) . Next, 
we use the inequaUty 



K+H-l 



k=K 

in order to conclude that 

K+H-l 





l+u 1 




V 271 


2n 


< 







K+H-l 



j^ Yl \S i9k) - Sx (gk) 



|2ra 



k=K 



1 

2n 



0(1) 



Y \S {9k) - SA9k)\^^ 



2n 



1+iL 

271 






„ , 271-1-1^ , ,„ 1 + v 1 + v 

O H^2^ (log log Kf-^ H^^ 



k=K 

Therefore, 

K+H-l 

Yj \Sx{9k)?'^~^~''\S{gk)- S^igk)]' 

k=K 



for 1 < z^ < 2n — 1. Hence, 

K+H-l K+H-l 

Y S{gkf^-' = Y SAgkf''-' + o(H{loglog Kr~' 



O i H (log log Kf'^ ) =0{H (log log K) 



\n-l 



k=K 



k=K 



0(H{loglogK) 



n-1 



D 



Corollary 6.2. Suppose 1 < ciK^ <H< C2K, where 6 > 1/2 and ci, C2 > 0. Then, (i) 



K+H-l 



- E 

H ^ 



k=K 



V27TS {tk + ^ak) 



Vlog log tk 



2n 



-^ ( 1 + ((log log i^y 

(27r) n! 



-1/2 



and (ii) 



K+H-l 



- y 



k=K 



'V2TTS{tk+M 

Vlog log tk 



2n-l 



O (log log i^) 



-1/2 



This Corollary implies Theorem 11.11 



7. COVARIANCE 

Lemma 7.1. Let primes Pi ^ P2 < V ^ cK. Assume 1 < H < cK, and let k' = k + x, 
where < x < K"^ , with e G [0, 1). Then, 

K+H-l / l/2i , \ 

Y exp{-iigklogpi- gk'logp2)) = O Ih^^j^^L + yy^K'/HogK\ . 
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Proof: By using Lemma [l!2l we can estimate: 



ln.^^ Jn.,..^^ '^VT ^^ , ^^J logy \ ^ J xXogy 



g" (t) \ogpi-g" (t + x) logp2 = — ; -^ (logpi - logp2)+o — ; To ]+o 



It follows that 



and 



\g" (t) logpi - g" {t + x) logpsl > c— — — 2 , 

yK (log K) 



\g" (t) logpi - g" {t + x) logp2| < cylogy- 



yK{\ogKf' 
The conclusion of the lemma follows by applying Theorem 13. II D 

Proof of Theorem 12.31 In order to compute EXJ X2 , we proceed as above in the 

calculation of eI'xJ^^) . 
Since by Proposition 14.11 

K+H-l 

E \S{gk)-S.A9k)\''' = 0{H), 

k=K 

therefore, it is essential to compute 

2Af-l 

jY X] Sx{gk)Sx{gk'), 

k=N 

where k' = k + a (log A^)^ . By using function ry, defined in ( fT5] ). we obtain: 

2N-1 2N-1 

Y^ Sx{gk) Sx{gk') = -——2 Y {ri{9k)r]{gk')-'n{gk)ri{gk')-ri{gk)ri{gk') + ri{gk)ri{gk')). 

k=N V^^'' k=N 

For the first term in this sum, we write 

2N-1 2N-1 

Y. ri {Qk) r] {gw) = Y. 7^^ Y. pVpV'"' ■ 
k=N Pi,P2<x^ yi'^f^ ^^^ 

Note that the sum 

2N~1 2N-1 

Y pV^'pT'^' = Y -M-^{gk\ogp,+g,^^^,^^^yAogp,)] 

k=N k=N 

is an exponential sum, and it can be estimated by using van der Corput's theorem by noticing 
that the second derivative of the function 

/3 



9 (s) log Pi +g \^s + a {log Nf j logp2 

is bounded by O ( N^^ (log N)^ j from below and by O ( N^^ (log N)~ log x ) from above. 
This implies that with an appropriate choice of x, 

2N-1 
T7 Y 'n {dk) 11 {gy) =0(1), 



N 

k=N 



-1 Y^2Ar-l. 



and similar for iV ^Y.k=N Vi9k)vi9k') 
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Therefore 



1 ^^"^ 1 

k=N 



2N~l 



^ Y Vi9k)vi9k'] 



k=N 



+ 0(1), 



where 



2N-1 



2N-1 



Y^ V{9k)r]{gk 



E -7^=^ E pi'^'pl'"'- 



k=N 



pi,P2<X^ 



VPIP2 



k=N 



If pi / p2, then by using Lemma 1770 we can estimate 



27V- 1 



Y - y— Y. Pi'^'pf = O (x^iVi/^ log at) = 0{N), 



P1,P2<X' 



VpiP2 



k=N 



provided that x = iV^ and k < 1/6. 
If pi = p2, then we have 

N z^ z^ p^ 

k=N p<x^ ^ 

If one sets x = N'^, then by using the definition of function 5, it is easy to see that for 

every k G [N, 2N - 1] , and k' = k + c (log iV)^ , we have 



gk>-gk = 27rc{logNf-^ + 0(^ilogNf-^ 

= (27rc/K(^-^^ (log a;)^-i + O ((log xf-^ 



where the implicit constant in the 0-term does not depend on k. 
Hence, by using the result of Lemma |7.2| below we find that 



2N-1 



h E E-p"^''"'"^ = (i-/')+iogi°g^+o(i) 



N ^^ ^^ p 

k=N p<x2 



It follows that 



2N-1 



^ E S, (gk) SAgk') = ^{1- 13) +loglogx + 0(1 



k=N 
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Next we note that 

2N-1 , 2N-1 



k=N k=N 



k=N 
2N-1 

+-^ XI -^^ (s-fc) {S {gk') - So: {gk')) 

k=N 
2N-1 

fc=Ar 
2N-1 



By the Schwarz inequahty, the last three terms can be estimated as (log log N) ' , and there- 
fore we have 

2N-1 



^ j; 5 (gk) S igk') = A (1 - /5)+ log log X + O ((log log 

k=N ^ 

This implies that 



D 

Lemma 7.2. Suppose that s (x) = c (log x)'^^ + O I (log x)^^ I , where /3 > 0. T/ien, 



y- \pis{x) = (1 _ /5) log log X + O (1) 

p<x 



Proof: Note that 

Y^ -pi<^) = Y^ n-i+i^^) [vr (n) - vr (n - 1)] , 

where vr (t) is the counting function for primes: vr (t) = |{p : p < i}| • Hence, by partial 
summation, 

Y V'^'^^ = vr (x) x-i+^^(^) - (-1 + is (x)) r TT (i) r^+^"(^)dt + O (1) • 

By using the estimate 

\ogt V(logt)V 

we obtain 

1 j-x is{x)\ogt 

Y-p'<^^ = {l-ts{x)) ^-^^dt + 0{l). 



18 



V. KARGIN 



Now, let US estimate the integral: 



I 



X is{x) logt 



dt 



/3 *logt 

If /3 > 1, then we can use the estimate 



and conclude that 



logx is{x)u 



-du = / — dv. 

logs ^ Js{x)log3 V 



f(a:)logx gTO 



-dv 



O 



X pis(x) logt 



tlogt 



dt = 



s{x) 



If/3G (0,1), then 



s(x) logx iv pi piv 

— dv = / — dv + 

s(x)log3 ^ Js(x)log3 ^ Jl 



s{x)logx iv 



-dv. 



(Note that s{x)logx — )• oo by assumption that /3 > 0.) The second integral is bounded 
uniformly in x. For the first integral we note that 



J six 



e'" -I 



{a;)log3 ^ 



and therefore 



J s\ 



-dv 



(x)log3 V 



So, if < /3 < 1, then we find that 

1 . 



dv = 0{l) 



1 



dv + 0{l) 

(a:)log3 '" 

-logs(x) + 0(1) 
l-/3)loglogx + 0(l). 



Y^ -pi<^) = (1 - /3) log log X + O (1) . 



p<x 



p 



If /? > 1, then we find that 



D 



y ipi^W = o (1) 

p<x 



8. Joint Moments 

Proof of Theorem |2.4t It is clearly enough to prove the corresponding result for ran- 
dom variables S {gki) and S{gk2) since X^ are the rescaled versions of these random vari- 
ables. In fact, as a consequence of the Selberg approximation result, it is enough to show that 
Sx idki ) and Sx (5^2 ) have the required moments. 
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Indeed, 

S {gkiT S {gk^f = {S^{gki) + S{gk^) - S^ {gk^)f [S^ {gt^) + S (gk^) - S^igk^))^ 
= Sx{gkiT Sx{gk2) 

+Oy^S^ {gk^y {S {gki) - S^ (5fci))""'' S^ {gk^Y (S {gk^) - S^ (fi'fca))''"* 

where the sum is over s and t such that 0<s<a, 0<t<6, and s + t < a + b. 
After we sum over ki and apply the Schwartz inequaUty twice, we find that 

2N-1 




2N-1 \ ' / , 2N-1 ^ ' 






ki=N I \ ki=N 



Hence, if variables 5 and Sx are scaled by (log log N) , the difference in their moments is 

1/9 

of order (log log N) ' . 

The result about moments of the scaled versions of Sx {gk, ) and Sx {gk2 ) follows from the 
result for random variables 

'^i •= /. 1 AT ^I [9kdN,L,)) , 

V log log A* 



where i = 1,2, and rj (t) is as defined in (fTSl ). 

Theorem 8.1. Let ai,a2,bi,b2 > 0. The joint moments of random variables 7]\ andu]^ , 

0.2 



mN{ai,a2,bi,b2) ■■='&\iii j 



WY^ ( JN)Y^ {JN)Y^ (JN) 



Vi \V2 \V2 



converge to the corresponding joint moments of complex Gaussian random variables rji and 
r]2, which have the following covariance structure: Ery? = Kfjf = E,riiri2 = ^riirJ2 = 0, 
^Virji = 1, Er?ir/2 = Efjirj2 = (1 - /3) , . 



Indeed, if this result holds, then the joint moments of (real) random variables 

converge to the corresponding joint moments of Gaussian random variables 5i and <S'2 , where 

E {S!) = E [Si) 



~^2\ _ TP f o2\ _ 1 

27r2' 
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and 



This implies the statement of Theorem [ 

Before attacking Theorem 18.11 let us recall the Wick Rule for the joint moments of Gauss- 
ian random variables, namely, 

7i-eP2({i,...,fc}) (r,s)eiT 

where the sum is over all pairings of indicies 1, . . . , A;. (In particular, if k is odd, then the sum 
is empty.) (See, for example. Theorem 22.3 in ||T2| or Appendix 1 on p. 13 in ifTSi ). 
If we apply this rule to random variables 7?i, %, then we find that 

m (ai, as, 61,^2) := E (^1)"^ (m)"' {mt m"' 

is zero unless ai + 61 = as + 62 . If ai + 61 = as + 62 , then 

m(ai, 02,61,62) = n(A;,ai,a2,6i,62)(l - /3)+ , (16) 

where n {k, ai, as, 61, 6s) is a number of ways to pair ai elements rji and 61 elements 77s with 
as elements ry^ and 6s elements r/s so that exactly k elements are connected with an element 
that has a different index. 

Also, we need a generalization of Lemma 13. 3 1 

Lemma 8.2. Let 

^logg/c -^logpfc I 
fc=l fc=l / 

/ a2+b2 ai+fei \ 

+g{x + u)\ ^ log^fe- ^ logPfc , 



,A;=a2 + l A;=ai+1 



w/iere g {x) is as defined in rfTTl) . {pi, . . . ,Pai+6i} 7^ {^1, • • • 902+62} and primes pi and qi 
are less than y < K"^ for all i and e < 1/n. Assume 1 < H < cK and u < a (log K)^. Let 
n = [(ai + as + 61 + 62) /2j . Then, 



K+H-l I „/Si 

Kl/2logi^ 






y ' logy 

k=K 

Proof: We can re-write the definition of h (x) as follows: 



h{x) = hi (x) + h2 (x) = g (x) 



(a2+f<2 ai+bi \ 

^ logqk- ^ logpfc I 
k=l k=l J 



02+62 ai+61 

+ (5(x + n) -5r(x)) j ^ loggfc- "^ logPfc 

,fc=a2+l fc=ai+l 



ON FLUCTUATIONS OF RIEMANN'S ZETA ZEROS 21 

The second derivative of the first term can be estimated as in Lemma[33J If x G [K, K + H] , 
then 

h'{ (x) G [A, k\] , 
where 

Cl 

A = 7T and k = C2v"' log y. 

y^K {log Kf 

For the second term, we note that 

{g{x + u)-g{x))" = g'"ie)u, 
where 6 £ [x,x + u] , and by using Lemma [X2l we find that 

h2 (x) = O r-^^^§f^ logy] = o{h'l (x) , 

provided that y < K^ with s < 1/n. It follows that h" (x) ~ h'[ (x), and the conclusion of 
the lemma follows by an application of Theorem |3T| as in Lemma 1331 D 
Proof of Theorem I8.lt By definition, we write 



E{,ry(r] (,f y^ Nf ^ 



Af(loglogiV)("i+"2+''^+''2)/^ 

',N-1 ( 

E E 



2N-\ 



kx=N \p<x2 ^^ 



E 



V 



bi 




where k2 = ki + 
t{p,p,q,q') : 



^ 2 Vp 

a (log A^) . If we expand the product of sums, we get for a general term 

1 (gi---ga2)*^'i(9a2+i...ga2+fe2r'=2 



^Pl...Pa^+h^qi---qa2+,^ {pi . ..pa^r"^ (Pai+1 • • • Pai+bi )'^'2 ' 

wherep:= {pi, . . . ,pa^) ,p' := (p^^+i, . . . ,^^^+^2), 9 := fe, •••,%!), andg' := (gb,+i, . . . j^b.+fej. 

By using Lemma [831 we find that after we sum this term over ki and divide it by 
N (loglog A^)^"^^'^^"'' ^~^ ^'' , we get a non-negligible contribution if and only if there is a 
pairing that puts every qi in a correspondence with a pj , so that qi = pj . In particular, it must 
be true that ai + 61 = 02 + 62 = n. Hence, the moment is asymptotically equivalent to 

. 2N-1 

N {loglog NT ^ ^ t{p,p',q,q% (17) 

y l^ >^ ^ ki=N p-p'=q-q' 

where p • p' denotes the product of primes in p and p', and similar- for q- q'.lfwe consider the 
sum over all {p,p') , in which at least one pi appears twice, then we can see that this sum can 
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be estimated as 



'2N-1 



o 



E Ei E^ =0(Ar(loglogiVr-), 



jyZ. ■ 1 .i / p 




which gives a negligible contribution to the moment. 

Otherwise, if every prime appears only once in {p,p') and if p ■ p' = q ■ q', then there is 
a unique pairing between elements of {p,p') and {q, q') . Let this pairing be called vr. That is, 
vr (i) = j means that pi = qj. 

The terms that satisfy pairing vr give the following contribution to the sum in (ITtT i: 

E^^ — ) [T.^—p — I +o{N{iogiogNr-') 

where 

"11 = |{«,j : 7r(i) = j, 1 < i < 02,1 < j < ai}| , 

^^12 = |{^,j : vr(i) = j, 1 < i < a2,ai < j < oi +6i}| , 

and so on. 

By Lemma IT2I this can be computed as 

((1 - /3)+)"^^+"^^ N (log log TV)" + 0{N (loglogiV)"-!). 

After summing over all pairings we find that the moment equals 

j;((l-/3)^)-^+"-+o((loglogiV)-i). 

7r 

Recall n {k, oi, 02, 61, 62) is a number of ways to pair ai elements r/i and 61 elements 772 
with a2 elements rji and 62 elements r/2 so that exactly fc elements are connected with an 
element that has a different index. That is, n {k, oi, 02, 61, 62) is the number of pairings tt for 
which ni2 + n2i = k. 

It follows that asymptotically, the moment tends to n {k, oi, 02, 61, 62) (1 — /3)^ , which is 
exactly the corresponding joint moment of the Gaussian variables that we obtained in formula 
(HHl. D 

9. Conclusion 

We have shown that the distribution of two zeta zero fluctuations approaches a two-variate 
Gaussian distribution. It is natural to expect that the distribution of any finite number of fluc- 
tuations approaches a multivariate Gaussian distribution with the covariance matrix KXiXj = 

(1 -/3ij),, where 



log\kj{N,uj)-ki{N,uj)\ 

^*-^ N^i^oo log log iV ' 



(3ij = lim 



and the limit is assumed to be positive and the same for all u. (The fluctuations are around the 
predicted locations ^^^(Af.aj)) and ki {N,io) are defined as in (|3]l with appropriate changes.) 



ON FLUCTUATIONS OF RIEMANN'S ZETA ZEROS 23 

The proof should proceed similarly to the proof of the result for two fluctuations; however, 
the notation becomes more cumbersome. 

The positive numbers Pij in the covariance matrix are not arbitrary but satisfy the ultra- 
metric inequality: 

/3ik < max{/3ij,/3jk} . 
More about processes with these finite-dimensional distributions can be found in Section 4 of 

m. 
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